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exposure, it will be possible to derive significant constraints on both core-collapse supernova models 
and the supernova rate in the near universe up to redshifts z < 2. 
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I. INTRODUCTION 

The cosmic background of neutrinos generated by core- 
collapse supernova explosions throughout the universe is 
known as supernova relic neutrinos (SRN) [llor more 
precisely as diffuse supernova neutrinos (DSN) [2| . These 
neutrinos are generally believed to provide a new source 
of information on the core-collapse supernova explosion 
mechanism and both on the supernova rate (SNR) and 
on the star formation rate (SFR) up to high redshifts of 

ma. 

However, the detection of DSN is demanding, as both 
their integral flux of ~ 10 2 i/cm _2 s -1 and their energy 
of E v < 50 MeV are low. In general, all v and v flavours 
are contained in the DSN. However, z7 e are the most 
likely to be detected as the inverse beta decay reaction, 
D e +p — > e + +ri, has the largest cross section in the DSN's 
energy region [l[ . Liquid scintillator detectors (LSD) pro- 
vide a large number density of free protons, a clear de- 
layed coincidence signal for this detection channel and 
high energy resolution at low energies. A large- volume 
LSD is therefore an ideal candidate for DSN detection. 
The best experimental limit on the DSN flux <£> Pc < 
1.2 cirrus- 1 for E Pe > 19.3 MeV (90% C.L.) has been 
achieved by the Super-Kamiokande experiment 0, Q. 
However, this limit is strongly determined by the back- 
ground events present in a pure water Cherenkov detec- 
tor (WCD). As a liquid-scintillator detector (LSD) allows 
far better background discrimination, it opens a nearly 
background-free energy window from ~ 10 to 25 MeV 
for DSN detection. The target masses of KamLAND 
and of near-future experiments like BOREXINO 0] and 
SNO+ Q are not sufficient to reach significant statis- 
tics [9|. A large-scale detector like LENA (Low Energy 
Neutrino Astronomy) [ljl EH with about 50 kt of liquid 
scintillator is required. 

We show in the present paper that LENA allows the 
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detection of ~ 10 DSN events per year in the almost 
background-free energy range from about 10 to 25 MeV. 
After an observation time of ~ 10 years a spectral anal- 
ysis of the DSN is possible, having implications on both 
core-collapse supernova (SN) models and the SNR up to 
redshifts z ~ 2. If no signal was detected, the new limits 
are significantly lower than all model predictions and im- 
prove the limit given by Super-Kamiokande by a factor 
of - 9. 

This paper mainly focusses on the observational window 
that is given by the remaining background of atmospheric 
and reactor v e that are indistinguishable from the ac- 
tual DSN signal. In Sections ITT1 and [Hl| a description of 
the planned detector and an overview of the predictions 
on the DSN flux and spectrum are given. In Sect. lIVI 
the v e detection channel is discussed. The results of the 
calculations on the site-dependent reactor neutrino back- 
ground rate and its spectral shape, considering especially 
the highly energetic part, are shown in Sect.[V] The cos- 
mogenic background, mainly 9 Li and fast neutrons pro- 
duced by muons, is shortly discussed in Sect. lVIl Using 
this information, the expected event rates for LENA in- 
side the energy window from — 10 to 25 MeV are given 
for different DSN models in Sect. lVIl! Furthermore, the 
spectral discrimination of different DSN models has been 
investigated by Monte Carlo simulations. Their results 
as well as a method to find a lower boundary for the SNR 
at z < 2 are also described in this section. Conclusions 
are given in Sect. I Villi 



II. DETECTOR CONCEPT 

The LENA concept foresees a cylindrically shaped de- 
tector of about 100 m in length and 30 m in diameter (see 
Fig. [l| . An inner part of 26 m diameter contains about 
5 x 10 m 3 of liquid scintillator, whereas an outer mantle 
with a thickness of 2 m filled with water is intended as 
a Cherenkov muon veto. The fiducial volume is defined 
within a radius of 12 m of the inner cylinder and encloses 
a volume of 4.4 x 10 4 m 3 . About 12 000 photomultipli- 
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FIG. 1: Sketch of the LENA detector 



ers (PMs, 50 cm diameter) installed on the inner surface 
provide a surface coverage of about 30%. The effective 
coverage could be further increased by mounting light- 
collecting concentrators to the PMs. 
At present, a preference exists for a liquid scintillator 
based on the organic solvent PXE (CieHig, p = 985 g/l) 
that was already tested in the BOREXINO prototype, 
the Counting Test Facility (CTF) at the Gran Sasso Na- 
tional Laboratory (LNGS) [l^ . Due to the large detec- 
tor radius, an attenuation length of ~ 10 m at 430 nm 
has to be achieved in order to obtain an adequate pho- 
toelectron (pe) yield in the PMs. It has been shown 
that this aim can be reached by p urification of the PXE 
in an Aluminum-column [l2l . Il3j . Furthermore, adding 
80 weight-percent of Dodecane (C12H26, p = 749 g/i) to 
the solvent has a positive impact on the transparency 
of the liquid and increases the number of free protons 
in the scintillator - and therefore the v e event rate - by 
almost 25 % [l]| . As Dodecane slightly lowers the light 
yield of the scintillator, both pure PXE and the described 
mixture provide similar photoelectron yields of at least 
lOOpe/MeV for an event in the center of the detector. As 
primary and secondary wavelength shifters (fluors) 6 g/i 
PPO and 20mg/£ bisMSB will be used [lj]. 
Besides detection of the DSN v e , LENA will be an ob- 
servatory for solar neutrinos [Io[, geoneutrinos 14 1, and 
the neutrinos emitted by a galactic core-collapse SN 
EE El EE EH- Moreover, the detector will allow to 
investigate the properties of neutrinos via beam experi- 
ments, and to search for proton decay [l7l [l8j. 
At the moment, the preferred detector sites are a mine, 
the Center of Underground Physics in Pyhasalmi (CUPP, 
Finland) [l9[, or the underwater plateau in the Mediter- 
ranean Sea used by the NESTOR Collaboration next to 
Pylos (Greece). Both sites provide an effective shielding 
of ~ 4000m.w.e. against cosmic radiation. We have cho- 
sen LENA at CUPP as our default scenario, as Pyhasalmi 
is far away from the middle-European reactors and is able 
to provide the required infrastructure for a large-scale ex- 
periment. In addition, local representatives have signaled 
big interest in the project. However, as the background 
due to power reactor /> e -flux is most important, locations 
in France, the US, next to the islands of Hawaii, and New 
Zealand have also been investigated for comparison (see 
Sect.©. 



III. THEORETICAL PREDICTIONS 

The energy-dependent DSN flux discussed in this 
pa per was taken from the publications by S. Ando et al. 
[l|,[2fj. In general, theoretical models for the expected 
DSN spectrum depend on two sources of input: 

Supernova neutrino spectrum. The spectral form 
of the DSN flux is strongly dependent on the supernova 
(SN) core-collapse model applied. Three different 
predictions have been made in [l[ with reference to the 
SN simulations performed by the Lawrence Livermore 
Group (LL) [2lj . by Thompson, Burrows and Pinto 
(TBP) and by Keil, Raffelt and Janka (KRJ) H|. 
The observed v spectra on earth will differ from those 
emitted by the proto-neutron star as the v and v are 
affected by the matter potential of the dying star [24], HE • 
As the neutrinos pass from the large matter potential 
of the central region to the surrounding potential-free 
vacuum, the change in the mixing parameters leads to 
a partial conversion of (p e into and vice versa. 



In 



In most SN models, the larger mean energy of ' 
at production hardens the resulting spectrum, 
addition, all models predict at least one resonant flavour 
conversion inside the SN envelope due to the mixing 
angle 612 that applies to the neutrino sector only. A 
second resonance might occur if the value of #13 is 
larger than ~ 1° and therefore enables an adiabatic 
conversion. This will affect either the v sector in case of a 
normal neutrino mass hierarchy, or the v for an inverted 
hierarchy. In the case of full adiabaticity the resonance 
will lead to an almost complete conversion of P P <-> v T as 
6»i3 is limited to values 13 < 12.5° [H, EE Hi- In thc 
following, we will refer to this resonant conversion of D e 
and v r as resonance case. All other combinations of the 
value of #13 and the mass hierarchy result in identical v e 
spectra. This situation we call no resonance case. 
If the mean energy of the v T is much larger than that of 
the v e as predicted by the LL model, the resonance case 
noticably shifts the spectrum of the v e to higher ener- 
gies. However, if one assumes only a small gap between 
those energies as in the KRJ model, the effect becomes 
negligible. For further explanation see fll. l24l [25|. 

Supernova rate. Additional information is re- 
quired on the redshift-dependent supernova rate Psn(z) 
(SNR) 1 . The SNR can either directly be measured 
[271 l28j or can be derived by combining the observed star 
formation rate p*(z) (SFR) with the initial mass function 
(IMF) of the forming stars [29(. As only massive stars 
with more than 8 solar masses (M Q ) will end in a core- 
collapse SN [lj|, their rather short lifecycles constitute a 



1 In the present paper, SNR is always referring to the rate of the 
core-collapse SN of type Ib,c and II. 
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FIG. 2: The z-dependent supernova rate calculated according 
to Ando et al. MJ (solid line) in comparison to the best fit 
(dashed line) [291 ] of direct SNR observations (diamonds) for 

z<iM,M. 



negligible aberration in comparison to cosmic time scales. 
Nevertheless, there is an uncertainty in the proportion 
of heavy SN-progenitor stars due to the uncertainty in 
the mass dependence of the IMF [29L l30l | at low masses. 
The standard Salpeter IMF (4>(m) oc m -2 35 ) returns 
the relation psn(z) — 0.0122 Mg 1 p*(z) [l|, whereas for 
modified IMF assumptions the conversion factor varies 
from —25 % to +8 % relative to the Salpeter result (30j . 



Resulting DSN spectra. In the present paper, we 
follow the model calculations by Ando [l[ in combining 
an observationally obtained SFR with the Salpeter IMF. 
We obtain a redshift-dependent SNR 



Psn(z) = 3.9 x 10 4 fsNh 7() - 
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fl m = 0.3 and J7a = 0.7 are the cosmological parame- 
ters of the matter and dark energy density, respectively. 
The overall normalization uncertainty in the SNR is pa- 
rameterised by fsN- The resulting values of the SNR 
f° r fsN — 1 are in reasonable agreement with direct SN 
observations for z < 0.9 [H [H, [29| , as shown in Fig.O 

In combination with the different SN models men- 
tioned above one obtains three models for the DSN spec- 
trum as depicted in Fig. [31 assuming /sjv = 1 and the 
absence of a matter resonance in the v sector. For the 
LL model, an additional DSN spectrum LLres is shown 
for the resonance case. 

Unlike the neutrinos that are produced by nearby SN 
the spectral portions corresponding to SN explosions fur- 
ther away will be increasingly redshifted by cosmic ex- 
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FIG. 3: DSN model spectra for v £ calculated by Ando et al. 
[2TJ ] using core-collapse SN simulations performed by the 
Lawrence Livermore group (LL) Thompson, Burrows 

and Pinto (TBP) [12] and Keil, Raffelt and Janka (KRJ) [1|. 
In addition, an LLres spectrum is shown in case of a resonant 
flavour conversion in the SN envelope. 
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FIG. 4: Contributions of different redshift regions to the DSN 
energy spectrum j20j. Neutrinos emitted by far away super- 
novae are redshifted by cosmic expansion. 



pansion [20j. For this reason, DSN from high-z regions 
will dominate the low-energetic part of the observed spec- 
trum. As shown in Fig.[U the flux of neutrinos from 
sources at z > 1 is substantially limited to energies 
E Pe < 10 MeV. 

As discussed above, the actual flux of DSN could vary 
widely from the predictions, as apart from the depen- 
dence on the used SN model there are uncertainties 
due to the assumed value of psn{z)- Direct measure- 
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ments reach no further than redshifts z = 0.9 and suf- 
fer from light extinction due to interstellar dust [29l |. 
The same is true for observations of star formation re- 
gions in the ultraviolet (UV) [U I321.J33T] and far infrared 
(FIR) [34], Hl[ band. According to [l|7 these uncertain- 
ties in the SFR translate to a factor Jsn between 0.7 
and 4.1 if one assumes a Salpeter IMF. However, com- 
bining the data of these observations with the cosmic 
gamma-ray background and the limits on the DSN flux 
by Super-Kamiokande, Strigari et al. have derived a con- 
cordance model of the star formation rate (CMSFR) with 
a favoured value of Jsn — 2.5 [36| . New measurements 
of the SFR further reduce the parameter space, as dis- 
cussed by Hopkins and Beacom [3(|. While these new 
predictions develop differently with redshift especially 
for z > 1, they are well contained inside the limits of 
0-7 < fsN < 4.1 given by Ando [J. 

Current UV and FIR astronomy and their future projects 
might be able to determine both the SNR and the SFR 
with high accuracy [4|. The present uncertainty of the 
SFR is only sa 30 — 50 % in the redshift region z < 1 
[30| . However, as the DSN detection is not suffering from 
dust extinction, neutrinos could test the validity of the 
assumed models for light extinction. Especially measure- 
ments in the energy region below lOMeV would provide 
valuable information on high redshift regions. However, 
in this energy regime reactor v e prove to be an undistin- 
guishable background and partially hide the DSN due to 
the high v e rates, as will be discussed in Sect.fVl 



IV. u e DETECTION CHANNEL 

As mentioned before, v R are the best choice for the de- 
tection of the DSN background. The cross section of the 
charged current interaction of a v & with a Hydrogen nu- 
cleus in the target, the inverse beta decay v e +P — > n+e + , 
is substantially larger than that of all other detection 
channels, a = 6.8 x 10~ 42 cm 2 at lOMeV [37[. More- 
over, the reaction provides a low energy threshold of 
1.8 MeV (corresponding to the mass difference between 
proton and neutron plus positron). 

Both decay particles can be detected in a liquid- 
scintillator detector, providing a coincidence that can 
be used for a very effective background reduction (see 
Sect.fV]) : The prompt signal due to the ionising processes 
and annihilation of the positron is followed by the signal 
of a 2.2 MeV gamma quantum. This is released after a 
delay of ~ 180 [is when the neutron is captured by a free 
proton in the scintillator, n + p —> d + 7. 
Due to the large mass of the neutron the positron will 
carry most of the energy that was deposited by the v e . 
In good approximation, the kinetic energy of the positron 
will be reduced by ~ 1.8 MeV compared to the incoming 
neutrino due to the Q-value of the reaction. However, 
the annihilation of the positron in the scintillator adds 
another 2m e c 2 of energy to the detected signal, leading 
to a total reduction of ~ 0.8 MeV. 



For the calculation of the expected event rate in LENA, 
the DSN flux must be convoluted with the energy- 
dependent cross section. The rate is then multiplied with 
the number of free protons in the fiducial volume. This 
number is 2.9 x 10 33 for the case of a scintillator mixture 
of 20 % PXE and 80 % Dodecane and a fiducial volume 
of 44 x 10 3 m 3 . 

We calculated the expected energy resolution in LENA 
by using the results of both laboratory measurements of 
the scintillator properties and Monte Carlo (MC) sim- 
ulations: The proposed scintillator mixture provides a 
light yield of ~ 8 x 10 3 photons per MeV, an attenuation 
length of ~ 10 m and a scattering length of ~ 30 m for 
the scintillation light (at 430 nm) [l3[ . Assuming a wall 
coverage of 30 % and a quantum efficiency of 20 % for the 
PMTs, a photoelectron (pe) yield of at least HOpe/MeV 
can be achieved for an event in the center of the detector 
[l7| . The resulting la energy resolution is 0.10/a/E (in 
MeV). 

For the three DSN models, we find in the no resonance 
case that LL [5l[ provides the largest event number: 
~ 6.8 fsN (see Eq. ([T])) detected D e in 1 year of measure- 
ment time in LENA, followed by the KRJ (~ 6.1 /sat) 
and by the TBP model (~ 4.7 Jsn)- This variation of 
rates originates from the different spectral forms of the 
models: As the cross section increases with energy, the 
LL model which predicts the largest flux above 10 MeV 
also provides the largest event rates. This is especially 
true for LL res (~ 7.7 /sn)- For all models, most of the 
9 e events are expected to be in the energy region between 
6 and 14 MeV. However, for these estimates background 
events have been neglected. 



V. v e BACKGROUND 

A large fraction of the v e events cannot be attributed 
to the DSN because of the background events due to 
v e generated both by air showers and nuclear power 
plants. These are intrinsically indistinguishable from the 
v e of the DSN and therefore independent of the type of 
detector used. 



A. Reactor ^ e 

For energies below ~ 10 MeV, the man-made back- 
ground due to nuclear reactors sets a threshold for 
DSN detection. The v e are generated by the (3~ decay 
of neutron-rich fission products of 235 U, 238 U, 239 Pu 
and 241 Pu inside the reactor. As the reactor v e flux is 
quadratically declining with distance, at least the closest 
reactors must be considered for a good estimate of the 
actual flux in the detector. 

Above the Q-value of 1.8 MeV, the spectral shape of 
the reactor neutrinos is best known up to energies of 
E Pe = 8 MeV, both from experimental data and from 
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TABLE I: Contributions of the fission products of uranium 
and plutonium to the fission processes and therefore to the v R 
flux emitted by a reactor (averaged over time) j39( ] . The last 
column shows the fission yields of the high-endpoint /3-emitter 
94 Br 0. 



fsotope 


Contribution to 
total fission rate 


94 Br 
fission yield 


235tj 


0.59 


f .66 x 10 _e 


238 jj 


0.04 


7.90 x f0~ 5 


239p u 


0.285 


2.7f x 10 -5 


241 Pu 


0.075 


L05 x 10" 6 



theoretical calculations that have been using the fission 
yields and (3 decay schemes of the isotopes in question. 
However, as the DSN flux is several orders of magnitude 
lower than that of the reactor neutrinos, it is necessary 
to take the high-energetic tail of their spectrum (up to 
E Pe ps 13MeV) into account when determining the lower 
detection threshold for DSN observation. 

Reactor Neutrino Spectrum. Especially for ener- 
gies above 8MeV, the exact spectrum of the reactor P e 
cannot be measured directly at a reactor because of the 
poor statistics. Instead, the spectrum is deduced from 
fission yields, [3 endpoint energies and decay schemes 
of the neutron-rich isotopes produced in a reactor. 
The experimental challenge is set by the extremely 
neutron-rich isotopes with high Q-values and lifetimes 
in the range of 10~ 2 seconds. 

Tengblad et al. have indirectly determined the reactor 
neutrino spectrum up to an energy of 12MeV by a 
measurement of the beta-decay spectra of the relevant 
fission products of 235 U, 238 U, 239 Pu These three 

elements are, on average, responsible for about 92% of 
the fission processes (see Tabled and therefore of the 
neutrino flux generated by the reactor [39| . In contrast 
to 238 U, 241 Pu contributes only a small portion of the 
fission products emitting high energetic neutrinos. For 
this reason, the spectral contribution of 241 Pu can be 
neglected to good approximation. 

However, there is at least one additional element known 
to be produced in fission processes which has an even 
higher (3 endpoint energy: 94 Br with a Q-value of 
13.3 MeV [13]. Due to its short lifetime of 70 ms and 
its low fission yields (as shown in Table|T|) its exact (3 
decay scheme is not known [4l| . Thus, in the present 
paper we have only been able to give an upper limit for 
its contribution to the spectrum, using its fission yield, 
Q-value, and the information that it decays in 70 % of 
all cases without emitting an additional neutron [4JJ. 

Site-dependent Reactor Neutrino Flux. Using 
the reactor v e spectrum just described, the reactor back- 
ground flux and the event rates in LENA were calculated 
at a number of different locations: At present Pyhasalmi 
(Finland) and Pylos (Greece) are the preferred detector 



sites. The Laboratoire Souterrain de Modane (LSM) at 
Frejus (France) was included as an example for a middle- 
European detector site. The US American sites Kimball- 
ton, Henderson, and Homestake Mine were used as they 
have uttered interest in a LENA-like liquid-scintillator 
detector and one of them will most likely be the home 
of the Deep Underground Science and Engineering Lab- 
oratory (DUSEL). Hawaii and New Zealand were chosen 
as they are far away from the nuclear power plants on 
the northern hemisphere and are from this point of view 
optimal detector sites for observing the DSN. 
When calculating the v e flux, in a first step the number 
of neutrinos emitted per second R„ c by a nuclear reactor 
can be derived: 

Re. =N v j iss -^e= (1.38 ±0.14) l0 20 P th[GW ] s"\ 

where N uJlss ~ 6 Q and E fiss = 205.3 ± 0.6 MeV [4^ 
are the average number of neutrinos and the mean en- 
ergy produced per fission, respectively, P t h is the thermal 
power of the reactor and e = 0.75±0.06 Q is the average 
fraction of time that a reactor is running (energy avail- 
ability factor e). The thermal power P t h can be found in 
the online databases of the International Atomic Energy 
Agency (IAEA) [3. 

In a next step, using the coordinates of each reactor and 
the detector site their distance d can be found. The flux 
without oscillations can be calculated dividing i? p<! by 
47rd 2 . 

Finally, for each distance to a particular power plant the 
effect of oscillations v e — > v^^ has to be taken into ac- 
count. This was achieved by convoluting the normal- 
ized spectrum F(E) with the energy-dependent oscilla- 
tion probability and weighing it with the expected inte- 
gral flux, 

<f>(E,d)dE = ■^F{E){l~sm 2 (29 12 )sm 2 {n£{E)/d))dE, 

(3) 

where £(E) is the oscillation length and #12 the solar mix- 
ing angle. This procedure was repeated for all reactors 
and the individual fluxes were summed for different de- 
tector locations. The necessary data concerning the coor- 
dinates was taken from the International Nuclear Safety 
Center (INSC) [45| and includes all major power plants 
running worldwide in the year 2005. Future changes as 
the launch of the Finnish nuclear plant TV03 as well as 
the possible shutdown of Swedish plants have not been 
taken into account but will not substantially change the 
v e fluxes. The integral v e fluxes and the corresponding 
event rates for LENA are summarized in TableQl] The 
calculated flux (<f> = 2.1 x 10 6 cm _2 s _1 ) in the Kamioka 
mine (Japan) is comparable to the one actually measured 
in KamLAND and is given for reference [4(| . In addition, 
Fig. [5] shows the energy spectra of the reactor v e at Fre- 
jus, Pyhasalmi and Hawaii, i.e. at a high, medium and 
low-flux site. 

For an estimate of the overall uncertainties, we have 
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FIG. 5: Spectra of reactor u e at Frejus, Pyhasalmi and 
Hawaii. Shaded regions correspond to experimental and 
model uncertainties. Above 12 MeV, we added an upper limit 
for the spectral contribution of 94 Br neutrinos. For compari- 
son, the DSN spectrum according to the KRJ model [23| with 
fsN = 1 (see Eq. fT}) is shown. 



TABLE II: Reactor v e fluxes and event rates for LENA at 
various detector locations. Uncertainties of both fluxes and 
event rates are ~ 14 %. 



Detector Location 


reactor u e flux 

(cm- 2 S - 1 ) 


reactor events 

(0.5 Mt yrs) 


Kamioka 


(J) 


2.1 x 10° 


2.5 x 10 s 


Frejus 


(F) 


1.6 x 10 6 


2.0 x 10 5 


Kimballton 


(US) 


6.4 x 10 5 


7.3 x 10 4 


Pyhasalmi 


(FIN) 


1.9 x 10 5 


2.1 x 10 4 


Pylos 


(GR) 


9.2 x 10 4 


11.5 x 10 3 


Homestake 


(US) 


7.5 x 10 4 


8.6 x 10 3 


Henderson 


(US) 


7.4 x 10 4 


8.4 x 10 3 


Hawaii 


(US) 


10.9 x 10 3 


12.4 x 10 2 


Wellington 


(NZ) 


5.4 x 10 3 


6.2 x 10 2 



taken into account possible deviations due to the frac- 
tion of annual runtime e as well as the temporal vari- 
ation of the abundancies of Uranium and Thorium iso- 
topes (~ 5%) [39]. In addition, the most recent values 
and uncertainties of the oscillation parameters Am 2 2 = 
8.2t° i x 10" 5 eV 2 and tan 2 6 12 = 0.40tg;$ 47] have 
been included in the calculations. In Fig.O these overall 
uncertainties are indicated as shaded regions. 



B. Atmospheric v e 

The flux of the atmospheric v e is increasing with their 
energy and starts to surpass the DSN signal at energies 
around 25 MeV. However, the total flux of atmospheric 



TABLE III: Dependence of the total atmospheric neutrino 
flux below 60 MeV on the detector location. The scaling factor 
Satm compares this flux to the one at the Kamioka site. 



Site 


latitude (N) 


Satm 


Hawaii 


1.5° 


0.8 


Kamioka, Pylos, Kimballton 


36.5°, 36.6°, 37.4° 


1 


Henderson, Wellington 


39.8°, -41.5° 


1.25 


Frejus, Homestake 


45.1°, 44.3° 


1.5 


Pyhasalmi 


63.7° 


2.0 



v e is dependent on the geographic (geomagnetic) lati- 
tude and will therefore depend on the detector site [48[. 
The individual spectra can be calculated using MC meth- 
ods. However, for easier comparison and in accordance 
with publications by Ando and Sato we have extrapo- 
lated the model spectra calculated by Gaisser et al. [48| 
(including the corrections by Barr et al. |49() to ener- 
gies below 60 MeV. We have estimated the atmospheric 
v e spectra for high geomagnetic latitudes as, for instance, 
for Pyhasalmi using the 3D simulations by Liu et al. [5fj| . 
As the energy spectrum of the atmospheric v e only mildly 
depends on location I48j , we left the spectrum unchanged 
and have only taken into account the dependence of the 
total flux on the detector site. Tablc lLUl shows our re- 
sults. The total flux of atmospheric v e at a particular 
site can be related to the flux at the Kamioka site by a 
scaling factor s a tm- 

In water Cherenkov detectors, atmospheric and 
provide an additional background source by creating "in- 
visible muons" [q"H with energies below the Cherenkov 
threshold. However, such muons do not pose a back- 
ground in a liquid-scintillator detector due to its different 
detection mechanism. 



VI. COSMOGENIC BACKGROUND 

Up to now, the discussion only included background 
events due to additional v e sources. However, muons 
that pass the detector or the surrounding medium (rock 
or water) also have to be considered. Most of the spal- 
lation products of these muons can be easily discrim- 
inated due to the signature of the inverse beta decay. 
Still, radionuclides like the /3-n-emitter 9 Li [52j or fast 
neutrons can mimic this e + — n coincidence. For an esti- 
mate of the rates, only events in the energy window from 
10 to 25 MeV have to be considered. All our calcula- 
tions were performed for LENA at Pyhasalmi assuming 
a depth of 3960 m.w.e. and a corresponding muon flux 
of (1.1±0.1) x 10- 4 /m 2 s H. 

A. 9 Li In-situ Production 

As the j3 endpoint of 9 Li is at 13.6 MeV, it will af- 
fect the DSN detection only in the lower energy region 
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of the observational window. For a rough estimate of 
the expected event rate, one can adopt the value derived 
for KamLAND: 0.6 events in 0.28 ktxyrs exposure for 
E > 9.5 MeV @. Scaling the mass to 44 ktxyrs fiducial 
volume in LENA, including a reduction of the integral 
muon flux by a factor of ~ 9 in Pyhasalmi (54[ and con- 
sidering the dependence of the production rate on the 
muon energy (oc E^ 75 ) [55| . the resulting rate is approx- 
imately 20 events per year. 

However, as the muon passes through the fiducial vol- 
ume, it can be clearly identified and a time as well as a 
volume cut can be applied. The 9 Li cannot travel far, 
as its half-life is T 1 / 2 = 0.18s. Therefore, excluding a 
cylindrical volume of 2 m around each muon's path for 1 
second (~ 5 x Ti 12) decreases the background sufficiently. 
As the fiducial volume is hit by a muon about every 5 s, 
one loses about ~ 0.2 % of exposure time. 



B. Fast Neutrons 

Whereas muon tracks in the inner detector or in the 
muon veto can be clearly identified, muons passing the 
surrounding rock generate a background of fast neutrons. 
Usually, such a neutron will be accompanied by a num- 
ber of charged particles that can be identified in the veto. 
However, there is a certain probability that the neutron 
will pass into the inner detector unnoticed. The neutron 
deposits its remaining kinetic energy in the scintillator 
and thereafter is captured by a proton, mimicking a real 
v e event. 

Using the estimates by Kudryavtsev et al. [fOJ for neu- 
tron production at 4000m.w.e. and assuming a mean ab- 
sorption length of 0.75 m [54j , one obtains a rate of about 
10 5 neutrons per year entering the muon veto. However, 
if one assumes a more realistic energy spectrum at the 
boundary between rock and cavern [54j | , simulations per- 
formed with Geant4 [56[ show that only ~ 7 neutrons 
per year will reach the fiducial volume and only 0.5 per 
year will generate a signal in the relevant energy region 
from 10 to 25 MeV. As mentioned before, this number 
will be further decreased by the detection of accompany- 
ing shower particles in the muon veto and is therefore a 
conservative value. In the further discussion we assume 
the neutron spectrum to be energy independent. 



VII. DETECTION POTENTIAL OF LENA 

This Section mainly deals with the dependence of the 
observational window and the signal-to-background ratio 
on the actual detector location. As Pyhasalmi is one of 
the preferred sites, a statistical analysis has been carried 
out for this place in order to investigate the sensitivity 
of the detector regarding the separation of different DSN 
model predictions and the possibility to give constraints 
on the SNR. The results of such a spectral analysis are 
described in part IVII CI 



TABLE IV: Lower and upper energy thresholds for DSN ob- 
servation and predicted signal and background events inside 
the energy window for LENA at different locations after 10 
years of measuring time and fsN = 1 (see Eq. (HJ). The range 
in the signal rates is determined by the LL and TBP model 
predictions (no resonance, see Sect . fill [I . The background in- 
cludes atmospheric and reactor z? e events and signals due to 
fast neutrons. 







Energy Window 


DSN S/B 


Detector Location 


(MeV) 


(0.5Mtyrs) 


Kamioka 


(J) 


11.1 - 28.1 


21-42/11 


Frejus 


(F) 


10.8 - 26.4 


22-41/12 


Kimballton 


(US) 


10.6 - 28.1 


23-44/11 


Pyhasalmi 


(FIN) 


9.7 - 25.1 


24-45/13 


Pylos 


(GR) 


9.4 - 28.1 


27-49/12 


Homestake 


(US) 


9.0 - 26.4 


28-49/13 


Henderson 


(US) 


8.9 - 27.2 


28-50/13 


Hawaii 


(US) 


8.4 - 29.0 


31-54/12 


Wellington 


(NZ) 


8.2 - 27.2 


31-53/12 




A. 


Energy Window 





Using the information on the two v e background 
sources in a LSD (see Sect. IV Al and lVBj) . it is possible to 
define an optimal energy window for DSN observation. 
As both the atmospheric and the reactor neutrino flux 
are dependent on the selected detector site, an individ- 
ual upper and lower energy limit for DSN observation has 
to be set for every location. In a spectral analysis it is 
desirable to include as much of the DSN signal as possi- 
ble, without increasing the background signal too much. 
We have chosen the lower (upper) energy thresholds at 
those energy values where the flux of the KRJ model 
{fsN = 1) begins (ends) to dominate in comparison to 
the total background flux (see Fig. [6] described below). 
The site-dependent detection thresholds are listed in Ta- 
ble llVl The best limits can be achieved at Hawaii with 
a window from about 8.4 to 29.0 MeV. At Pyhasalmi, it 
will be from around 9.7 to 25.1 MeV. 



B. Event Rates 

The detection thresholds as well as the event rates 
vary with the location that is chosen for LENA. As TBP 
provides the lowest event rates and LL the highest, their 
rates are given in Table HVl for the different detector 
sites. In addition, the number of background events 
inside the energy window is shown. 

At Pyhasalmi, we have obtained between 24fsN (TBP) 
to 45 fsN (LL) events in 10 years of measurement within 
the energy window (no resonance, see Sect. lIIl[ ). In case 
of LLres, 53 fsN events would be detected. For a value 
of fsN = 2.5 [3a], one expects therefore ~ 100 events in 
this time period. The background due to reactor and 
atmospheric neutrinos within the same time would give 
~ 8 events, and up to 5 events from fast neutrons have 
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FIG. 6: Event rates of reactor, atmospheric and DSN v e (LL, 
KRJ, TBP, see Fig.[3]l as expected for LENA in Pyhasalmi 
after ten years of measurement and fsN = 1. The shaded re- 
gion represents the uncertainties of the DSN rates due to fsN- 
The energy window is chosen such that the flux of the KRJ 
model exceeds the background flux. The Super-Kamiokande 
limit is also indicated. 



TABLE V: Using the calculated background levels in two en- 
ergy bins, an upper limit on the DSN flux can be achieved. 
Values are given for 10 years of measuring and 100% detec- 
tion efficiency both for signal and background. For compar- 
ison, lowest model predictions (fsN = 0.7 in the SNR) and 
the current Super-Kamiokande limit are also given. 



Energy Range (MeV) 


10.5-19.3 


19.3-25 




atmosperic v e events 


2.2 


3.6 


in lOyrs 


reactor v e events 


0.2 







fast neutron events 


2.7 


2.0 




total 


5.1 


5.6 


in lOyrs 


flux limit (90% C.L.) 


0.3 


0.13 


cm^s -1 


lowest model prediction 


1.4 


0.16 




Super-Kamiokande limit 




1.2 





tions for both energy regions. 
(0.13 cm^s- 1 ) obtained by LENA are 20 % (80 %) of the 
flux for fsN = 0.7 below (above) 19 MeV. In addition, the 
limit from Super-Kamiokande is shown. LENA could im- 
prove the present limit by a factor of 9 within 10 years. 



C. Spectral Analysis 



to be added. Taking into account the overall uncertainty 
in the supernova rate for z = (factor fsN = 0.7 — 4.1, 
see Eq. (JTJ) the predicted DSN event range is widened 
further to 17 — 220 events. However, the direct upper 
limit on the DSN flux of Super-Kamiokande Q corre- 
sponds to 185 events in LENA. In any case, a significant 
detection signal of the DSN appears to be certain. The 
spectra of both signal and background events are shown 
in Fig. [6] for LENA located at Pyhasalmi. 

It should be emphasized that in Pyhasalmi about 
25 % of the registered DSN in the energy window are 
originating from redshifts z > 1. Their flux is expected 
to influence the spectral form at low energies visibly. 

Flux limits. Assuming there is no signal of the DSN, 
we calculated the flux limit that LENA would be able to 
give for two energy regions: One for the most sensitive 
region from 10.5 to 19.3 MeV where the background lev- 
els are lowest. A second one in the regime from 19.3 MeV 
to 25 MeV that corresponds to the limit achieved by the 
Super-Kamiokande detector [f|. Assuming a 100% effi- 
ciency for the detection of all v e and other background 
events, one can derive the statistical uncertainty corre- 
sponding to 90% C.L.. If one assumes the spectral form 
of the KRJ model, one can match the number of DSN 
events to the uncertainty of the background by scaling 
the DSN flux by means of the factor fsN- The derived 
limits on the flux are well below the current model pre- 
dictions. TableM shows background event rates, flux lim- 
its after 10 years of measuring and lowest model predic- 



According to the current DSN model predictions, 
LENA will provide a sufficient event rate for probing 
both parameters that enter into the actual DSN spec- 
trum: the redshift-dependent supernova rate (SNR) and 
the core-collapse SN v e spectrum. 

Limits on SNR. The event numbers in the energy 
regime from 10 to 14 MeV are quite similar for all DSN 
models, yielding between 1.48 fsN to 1.97 fsN events per 
year, as can be seen in TablefVT] Using the actually de- 
tected event rate in this energy bin, to derive limits on 
the scaling- factor fsN in the SNR the SN neutrino spec- 
trum does not have to be known. Fig. [7] shows LENA's 
exclusion potential if one assumes fsN = 2.5 to be true. 
The depicted curves show the significance level for re- 
jection of other fsN values dependent on the exposure. 
After 10 years of measurement, Jsn < 1-3 could be ex- 
cluded at a 2cr level. 

In this way, LENA will be able to independently ver- 
ify the value of fsN derived from optical observations of 
the SNR, and cross-check the validity of the applied dust 
corrections. At present, SNR observations reach out to 
z ps 0.9 [13, [H, . With the lower energy threshold of 
9.7 MeV for LENA in Pyhasalmi, there will be also a non- 
neglegible amount of DSN events due to SN at 1 < z < 2 
that can be used to derive limits on the SNR for this 
redshift region. Combined with assumptions about the 
IMF [3(|, DSN observation will also be able to provide 
constraints on the SFR. 

Discrimination of DSN models. The progress in 
optical observations is likely to provide a solid prediction 
for both the SFR and SNR up to z w 2. Using this infor- 
mation, an analysis of the DSN event spectrum in LENA 
could be used to constrain the parameter range for dif- 
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TABLE VI: Rates for the DSN models used in the MC- 
calculations and x 2 -tests performed in order to analyse the 
potential for model discrimination. Values are given for 
50ktyrs of exposure (one year of measuring time in LENA) 
at Pyhasalmi and fsN = 2.5. 





Rates (Nbi) 


Rates (N B 2) 


DSN model 


9.7-14.5 McV 


14.5-25.1 MeV 


LL (no res.) 


4.93 


6.20 


TBP 


3.70 


2.65 


KRJ 


4.88 


4.73 


LL (res.) 


4.50 


8.70 


background 


0.60 


0.73 



ferent SN models or even to discriminate between them. 
For testing LENA's potential, MC simulations of typi- 
cal DSN event spectra have been carried out. In these 
tests, one DSN model and the currently favoured value 
°f Jsn = 2.5 have been assumed to be true, and in ac- 
cordance with these assumptions 10 4 MC spectra have 
been simulated. As the event numbers are low, only 
two energy bins have been chosen: the first one reach- 
ing from 9.7 to 14.5 MeV and the second one from 14.5 
to 25.1 MeV. In Table lVTl the expected event numbers 
per year for the three DSN models and the background 
rates are shown, including the resonance case for the LL 
model (see Sect. lIIII) . As only two energy bins have been 
chosen, the excellent energy resolution of the detector of 
better than 3 % for Ep e > 10 MeV and the reaction kine- 
matics have almost no effect on the analysis. Both have 
therefore been neglected. 

Via a x 2 -analysis, an exclusion probability can be given 
that a MC spectrum created according to a given model 
is wrongly assigned to a different combination of event 
numbers and therefore to an alternative SN model. Fig. [5] 
shows exclusion plots for each of the three models assum- 
ing an exposure time of 10 years. 2a and la exclusion 
regions are located outside the depicted curves. For com- 
parison, the predicted event numbers for all models are 
shown. 

The evaluation shows that the separation potential for 
the models LL and TBP is best, as their spectral slopes 
differ most. For Jsn—2.5, LENA is able to discriminate 
between these two models with a significance of 2.6a af- 
ter 10 years of exposure. However, a separation at 2a 
level between spectra with very similar slopes as KRJ 
and LL would require long exposure times of > 30 years 
unless fsN is very close to the upper limit (fsN < 4.1). 
For the case of the LL model, the exclusion probability 
for the other models is shown in Table lVlTI as a function 
of the exposure time. Note that an exclusion at > 2a of 
a resonant flavour conversion in the SN envelope (corre- 
sponding to LLres) would be possible after 15 years of 
measuring. 




0) 

TO 2 

o 




exposure time (yrs) 

FIG. 7: Potential to derive bounds on the scale factor fsN in 
the SNR. If one assumes fsN ~ 2.5 according to the CMSFR, 
the curves shown in the plot depict the significance level at 
which other values of fsN can be rejected after a certain ex- 
posure time. 



100 




20 40 60 80 100 
event number, 15-25MeV 

FIG. 8: Exclusion plot for the assignment of a simulated 
event spectrum in LENA to a wrong DSN model. A value of 
fsN — 2.5 and 10 years of exposure are assumed. MC spec- 
tra created according to one of the models are compared to 
all possible combinations of event numbers in the two energy 
bins (see Tab. lVI[l . Regions of more than la and 2a exclusion 
probability for a wrong assignment are located outside the 
depicted lines. Predictions assuming the no resonance (dots) 
and the resonance case (diamond) are shown. LL and TBP 
model can be discerned at a significance level of more than 
2.6a. 
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TABLE VII: Exclusion probability for a wrong model assign- 
ment of a simulated DSN event spectrum in LENA. The LL 
model [no resonance case) is assumed to be true, significance 
levels of the rejection of other models are shown as a function 
of the exposure time and fsN- 



Parameters Significance of exclusion (a) 



fsN 


exposure (yrs) 


TBP 


KRJ 


LL res 


2.5 


5 


1.8 


1.1 


1.5 




10 


2.6 


1.3 


1.8 




15 


3.0 


1.5 


2.1 




20 


3.3 


1.6 


2.4 




2:> 


3.7 


1.8 


2.6 




30 


>4 


1.9 


2.9 


4.1 


30 


>4 


2.3 


3.6 



VIII. CONCLUSIONS 

It has been shown in this paper, that in a large 
liquid-scintillator detector (LSD) like LENA diffuse su- 
pernova neutrino (DSN) events can be detected in an al- 
most background-free energy window between ~ 10 and 
25MeV. High background suppression is achieved due 
to the fact that the neutron produced in the detection 
reaction - the inverse beta decay - can be detected in a 
LSD. Above ~ 10 MeV the high-energetic tail of the man- 
made reactor D e is becoming negligible in comparison to 
the DSN flux. Above ~ 25 MeV the atmospheric v e flux 
starts to dominate. 

If placed in Pyhasalmi, the lower threshold for LENA will 
be close to 9.7 MeV. For this location and for the most 
likely value of fsN — 2.5 (see Eq. ([T])) in the supernova 
rate (SNR), about 6 to 13 events per year will be con- 
tained in the energy window and can be detected. This 
would be the first detection of the DSN background. It 
should also be noted that about 25 % of the detected DSN 
will be originating from a red-shift region 1 < z < 2. If no 
signal was detected, a new limit of 0.13cm _2 s _1 on the 
flux above 19.3 MeV could be achieved within 10 years 
that would surpass the one of the Super-Kamiokande de- 
tector by a factor ~ 9. In the lowest background re- 
gion between 10.5 and 19.3 MeV the limit of 0.3cm" 2 s~ 1 
would be 20 % of the lowest current model predictions. 
Apart from mere detection, LENA will be able to distin- 
guish between different DSN models and give constraints 
on the form of the neutrino spectrum emitted by a core- 
collapse supernova. This can be reached via an analysis 
of the DSN's spectral slope. The significance of the re- 
sults will be highly dependent on the exposure time and 
the supernova rate in the near universe. For a known 
SNR with fsN = 2.5, the discrimination between the 
discussed LL and TBP models for the DSN will be pos- 
sible at a 2.6cr level after 10 years of measuring time. 
Distinguishing between DSN models with more similar 
spectral slopes, however, would require higher statistics. 



In addition, by an analysis of the flux in the energy re- 
gion from 10 to 14 MeV the SFR for z < 2 could be 
constrained at high significance levels. The current lower 
bound of fsN — 0.7 [lj could be increased to fsN = 1-3 
at a 2a level within 10 years of measuring time. 
LENA will therefore not only be able to detect the DSN; 
it will also be able to make valuable contributions to 
both, core-collapse SN models and the redshift depen- 
dent supernova rate. 
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